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Treating a many-body Fermi system in terms of a single 
particle in a deforming mean-field, we relate adiabatic geo- 
metric phase to susceptibility for the non-cyclic case, and to 
its derivative for the cyclic case. Employing the semiclassical 
expression of the susceptibility, the expression for geometric 
phase for chaotic quantum system immediately follows. Ex- 
ploiting the well-known association of the absorptive part of 
susceptibility with dissipation, our relations may provide a 
quantum mechanical origin of the damping of collective exci- 
tations in Fermi systems. 
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Chaotic adiabatic systems |]J serve as useful models for 
complex systems in the mean-field approximation. These 
idealized models have been employed to understand very 
interesting phenomena in nuclear Q and plasma physics 
H- In particular, an important problem of many-body 
theory is to relate the slow collective excitations to faster 
single-particle motions. Adiabatic approximation leads 
us to linear response theory on one hand where dynamical 
susceptibility (or polarization propagator) jl| is central, 
and geometric phases on the other. A relation between 
these seemingly disparate quantities is being sought for 
in this Letter. 

Phase factors of geometric origin were discovered for 
adiabatic quantum systems ^,[1 and have been gener- 
alised to non-adibatic situations [p| . The geometry of the 
Hilbert space plays a key role in understanding geometric 
phases ||, a recent illustration being for the Josephson 
junction Adiabatic geometric phase led to an under- 
standing of several phenomena in physics |l0| like frac- 
tional statistics in two-dimensional statistical mechanics, 
integer quantum Hall effect, anomalies in field theory, 
the Magnus force in the context of superfluidity and 
so on. Recently, it has been shown that the viscosity of 
quantum Hall fluid in two dimensions at zero tempera- 
ture is related to adiabatic curvature (whose flux gives 
the phase) An interesting aspect of many-body 

physics related to these advancements is in the Born- 
Oppenhcimer approximation where reaction forces are 
shown to be of geometric origin |13|| . Owing to an in- 
timate tangle of collective and single-particle effects, it 
becomes useful to treat the total system as one where 
a slow subsystem is coupled to a faster one. In this 
setting, first-order corrections to the Born-Oppenheimer 
approximation leads to geometric magnetism and deter- 
ministic friction in a fully classical treatment, whereas a 
half-classical treatment gives rise to geometric magnetism 
only. 

When a particle (a nucleon, say) moves inside an enclo- 
sure whose boundary is adiabatically vibrating in time, 
the wavefunction can acquire a geometric phase over a 
cycle of vibration. This model is an idealization of a 
single particle in a mean-field represented by the enclo- 
sure, and has been an established paradigm for numerous 
studies in the past |l| . We would particularly like to em- 
phasize here the importance of such studies in enhancing 
the understanding of damping of collective excitations in 
nuclear physics. Concepts like time correlations and sus- 
ceptibility are fundamental to any discussion of quantum 
statistical mechanics as they lead to an understanding of 
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the response of many-body system m and the transport 
coefficients [ fl4[ . It is well-known that the imaginary part 
of susceptibility is related to dissipation, thus it may be 
quite interesting to explore its possible relationship with 
a quantity like geometric phase. To work out such a rela- 
tion, we begin with a discussion of geometric phase and 
identify the appropriate quantities related to susceptibil- 
ity. 

Let us first consider a Hamiltonian parametrised by 
R which describes a single particle in an effective mean- 
field where the shape of the field is vibrating adiabat- 
ically in time. It is well-known |j| that when the pa- 
rameters evolve along a cyclic path, C, the instantaneous 
eigenfunction of the system, |n(R)) corresponding to the 
eigenvalue i?„(R), acquires a geometric phase given by 

7n (C) = I z(n(R|V R n(R)).dR = ~ f V n .dS, (1) 

where S is the surface enclosed by C in the parameter 
space, and V„ is the "field-strength" (adiabatic curva- 
ture) given by a familiar expression involving a sum of 
weighted wedge product between two appropriate matrix 
elements : 



V, 



(n|V R JT|m) A (n\V R H\m) 



(E n — E m ) 2 



A suitable form of V„ for the sequel is JTf 



V„ = — lim 

2h e-*o 



d<e- £t i(n|[(V Rj ff) t ,A(V R i?)]|n) (3) 



where (V R i/)t denotes the Heisenberg-evolved operator. 
Note that, the state |n(R)) appearing in (3) corresponds 
to a single-particle eigenket in an effective mean-field. 
This state is clearly related to the original many-body 
Fermi system for which the imaginary part of the dy- 
namical susceptibility is Q 



X (t) = -<$oP(*),£(0)Po> 



(4) 



where | <&o } is the pure ground state of the many-particle 
system with Fermi energy, Ep. If one-body operators, H, 
A, and B, are used to construct many-body operators by 
a direct sum so as to get H, A, and B, respectively, and 
H\$l) = then we @ have 



{$o\A\$i) = (m(R)|i|n(R)>. 



(5) 



On reducing the many-body system at T = OK, (where 
the Fermi-Dirac distribution is a Heaviside step func- 
tion), to one-body system, we can express Jig] 



dte iwt tr 5{E F - H)[A(t), B(0)]. (6) 



This can also be written semiclassically as 



2^h 



% |det(m;-I)|V2 ' 
. J dte lut <j> <1tA w (t)B w (t +t)+O(h ) (7) 

where the subscript 'W refers to the Weyl symbol of 
the operator, / denotes the degrees of freedom and L c i 
is the Liouvillian operator. The last term corresponds 
to the periodic orbit expansion where S p , v p , and m p 
correspond to action, Maslov index, and the monodromy 
matrix for the periodic orbit, p, and r denotes the repeti- 
tions. We have used the Gutzwiller trace formula |[L7| for 
the case where the single-particle motion is chaotic. The 
above semiclassical expression is valid for Hamiltonian 
operators which are quadratic in momentum, p, with an 
additive position-dependent term. The expression (6), 
however, is general. For the case where the system has 
partially broken symmetry, the results have been recently 
generalized fl§}| . 

The label V in (3) corresponds to single-particle states 
and is related to |$o) because the many-body matrix el- 
ements can be written in terms of one-body matrix el- 
ements for the case when all the constituents are taken 
as non-interacting. In many-body physics, this gives the 
zero-order response whereupon the interaction can be in- 
cluded in a Vlasov description in an iterative way p9| . 
For relating the response function to the geometric phase, 
the operators A and B in our discussion are to be iden- 
tified with V X H and V Y U for R = (X, Y, Z). 

The matrix element in (3) can be written as a many- 
body matrix element using (5) by composing A and 
B so that we get the operator, C(t) — [A(t),B(0)] — 
[B(t),A(0)], which s related to a difference x" AB(t) — 
X" BA(t) = x v c(t)- Thus, we can re- write V„ as 

V„ = — lim/ dte^ t (* |C|$o> 



dt t x "c(t) = 



du 



»=o- 



(8) 



We now arrive at our first relation for the case of cyclic 
evolution : 



7n(C) 



rfS 



9x"c^;R) 



du> 



(9) 



Since x"c(^) is related to energy dissipation, this rela- 
tion connects geometric phase to dissipation. 
Exploiting (7), we get 
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5 £ib( w ; r ) 



[A w (e" 



Cdt B w ) 



If f d f xd f p c/ 

— dt £6(E F - H) 

2 J J (2nh)f 1 ' 

B w (e- c <**A w )] 



- — lim 6{Ep -H)^ | det(m . _ I) |i/ 2 



die 11 



dr [Aw(r)Bwr(T + 1) - B W (Y).V(Y + *)] ( 10 ) 



which entails the semiclassical expression for adiabatic 
geometric phase for chaotic systems, using (9). The first 
term comes from the classical time correlation function 
which is expected to decay exponentially as the dynamics 
is chaotic. This decay is governed by the largest Liapunov 
exponent of the classical system. The second term is the 
semiclassical correction in terms of periodic orbits which 
can be termed as a periodic orbit two-form. It is impor- 
tant to note that (9) is valid for general Hamiltonians 
whereas (10) is restricted to the operators like p 2 + V(q) 
only. 

The time-dependent deformation of the mean-field can, 
in general, be non-cyclic, particularly due to the fact that 
several harmonics (possibly incommensurate) may be in- 
volved. In this general case, the above expression does 
not hold and we now proceed to quantify the general re- 
lation. 

When a quantal system does not follow a cyclic evolu- 
tion the geometric phase not only depends on the path 
of the evolution curve but also on the end points. It has 
been shown in pfj that the wavefunction of the system 
can acquire a geometric phase during its time evolution, 
irrespective of adiabatic and cyclic condition, which is 
given by 



7 (r) = / i(x(t)\x(t))dt, 



(11) 



where |x(t)) is a modified state vector defined from the 
actual state vector \^(t)) of the system as \x(t)) — 
fflfflj-W) and an overdot denotes the timc- 



|(*(t)|*(0)> 

derivative. From this general expression, the adiabatic, 
open-path geometric phase can be written in the follow- 
ing way, 

7«(r) = jT*<Xn(R)|Vx„(R)>.dR = jTn„(B.).dR, (12) 

where f2„(R) is a generalised adiabatic vector potential 
(connection one-form) whose line integral gives the geo- 
metric phase. This non-cyclic adiabatic geometric phase 
can be expressed as a line integral of the difference of two 
vector potentials pl| , 



7„(T) = J(A n (R) - P„(R)).dR, 



(13) 



where A„(R) = i(n(R)|Vn(R)) is the usual Berry po- 
tential whose curl gives the adiabatic curvature V„ and 



P„(R) = -S( ^"n(R(o))^(R)^ ) is an extra potential that 
takes care of the contributions from the end points of the 
evolution path. The non-cyclic geometric phase is gauge 
invariant because under a local gauge transformation, A„ 
and P„ transform in the same way. Also in the special 
case of cyclic evolutions of parameters, (12) goes over to 
the well known expression, (1). 

For our purpose, we simplify the generalised vector po- 
tential, fi n (R), as 



Q„(R)=A„(R)-P n (R) 

(n(R(0))|m(R)> (m(R)\VH\n(R)) 



i(R(0))|n(R)> 



Using an integral representation of (E-, 
can write f2„ (R) as follows 



(E n — E m ) 

E m )~ l , 



(14) 



we 



- lim / 



dte- et $ 



te 



-i(E n -E m )t/H 



<n(R)|P n (R(0))P m (R)Vff|n(R)) 



|(n(R(0))|n(R)>|* 



(15) 



where P n (R(0)) = |n(R(0)))(n(R(0))| and P m (R) = 
|m(R))(m(R)| are the projection operators correspond- 
ing to nth and mth eigenstates. Now define a quantum 
correlation function between the operator, A, and the 
Hciscnbcrg evolved operator, B t , through 

QAB(t) = -(n\(AB t + B t A)\n) - (n\A\n)(n\B t \n) (16) 

and identify the Hermitian operators A = P n (R(0)) and 
B = VH. With this, the generalised vector potential can 
be written in terms of a quantum correlation function: 



On(R) = T lim 



dte 



1AB 



(*) 



|(n(R(0))|n(R))|s 



(17) 



Let us note that QAs(t) defines fluctuations in the sym- 
metrised correlation of two operators. Recalling the re- 
lationship between the Fourier transform of Qab (t) and 
the susceptibility x" (u) at zero temperature, we find 



>AB 



(w) = W{^\ 



(18) 



an instance of the fluctuation-dissipation theorem of the 
second kind p3] | . Following the physical arguments used 
to express single particle expectation values in terms of 
statistical quantities as in the cyclic case, we get 
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X "(0;R) = -lim 

= fi«(R), 

»(*) 



dte~ et Q' 



AB 



(*) 



(19) 



where Q' AB (t) — [^(R^ff^R)^ is the scaled quantum 
correlation function. Here, the scale factor is the survival 
probability P(t) = (ra(R(0))|n(R)) | 2 which is related to 
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the two-time correlation function at zero temperature for 



a system in canonical equilibrium 24 1 . This result seems 
to be a fundamental one as the generalised adiabatic vec- 
tor potential is exactly equal to the imaginary part of the 
susceptibility at zero frequency. Therefore, the adiabatic 
geometric phase for a non-cyclic variation of external pa- 
rameters along a path T is given by 



7n(r) 



^F(0;R).dR, 



(20) 



which relates a bulk property of the many-body 
Fermionic, chaotic quantum system to a physical quan- 
tity assigned to a single particle. 

In (9) and (20), note that the frequency-dependent re- 
sponse functions are different. However, they are related 
as follows : for a cyclic evolution along a path, C, apply- 
ing Stokes' theorem, (20) reduces to a surface integral of 
curlx"(0;R) over S. Using (9), it follows that 



V R xx"(0;R) 



dx"c(^;R) 



duo 



(21) 



The above results have clear relevance to finite Fermi 
systems like metallic clusters and nuclei. The historical 
collective picture of nucleus was Bohr's liquid drop which 
has been successfully used ever since. Due to the nu- 
merical simulations [Q and the results (eqs.(9) and (20)) 
obtained above, it is clear that the nuclear fluid has, in 
fact, characteristics of a gel which is an elastic or a vis- 
coelastic liquid. The connecting relation (21) between 
cyclic and non-cyclic case is necessary as a consistency 
requirement. The applicability of our results to nuclei 
is due to the compelling evidence, originating from shell 
structure of binding energies, testifying that the mean 
free path of nucleons is larger than the nuclear dimen- 
sions Hence it is natural to look for the source of 
viscosity in the nuclear fluid in terms of single-particle 
features where the effective mean-field leads to chaotic 
dynamics. 

It is well known that viscosity or friction in a quan- 
tal context can appear in thermodynamic limit. Since 
we have consistently discussed many-body Fermi systems 
in equilibrium, the seemingly contradictory conditions of 
finiteness of the system and "continuity" of the energy 
spectrum are met with. The finite-size of the system ex- 
plicitly manifests itself in terms of the sum over periodic 
orbits. Thus, yet another way to look at our relations 
(9), (20) is that these relations are off-springs of finite- 
size and spectral continuity- a stylized version of thermo- 
dynamic limit where viscosity and geometric phase can 
be discussed together. 

The notion of viscoelasticity in the context of nuclear 
physics has been pointed out in Ref. [1]. If we accept 
that the model for nucleus is a gel model and apply lin- 
ear viscoelastic theory |26 where stress relaxation func- 
tion is connected to the time-correlation function, we can 



see that the generalised vector potential (connection one- 
form) is directly related to viscosity. This is due to the 
fact that the time integral of the relaxation function is 
identified with viscosity in linear viscoelastic theory. 

Every time some relationships with possible fundamen- 
tal relevance are given, the question about their gener- 
ality arises. The relation (9) and (20), so-to-say (j n -x) 
relations, are restricted by the limits T — > 0. For any 
non-zero T, the semiclassical result for susceptibility is 
already known [|l6|. We believe that the (7n~x) rela- 
tion will be given as above where x" will be evaluated at 
some non-zero value of to consistent with the tempera- 
ture. Hence, as long as the dynamics is chaotic, these re- 
lations hold. This takes us to a deep interplay of quantal 
and statistical features embodied in the subject of quan- 
tum chaos. It is interesting to mention here that recent 
works have led to a new treatment of quantum statistical 
mechanics in terms of Berry conjecture p7j , f28f . Further, 
it has been shown ^] that thermal fluctuations for large 
number of degrees of freedom coincide with quantum fluc- 
tuations for chaotic quantum systems in the semiclassi- 
cal limit when the freedoms are smaller in number. One 
can immediately rationalise above discussion in the light 
of these contemporary works. In the above discussion, 
when we pick up a single-particle state, it ought to be 
one that can be written as a random superposition of 
plane waves in the sense of [ p7| . This is possible for a 
system with chaotic classical dynamics. It is these states 
which can be combined into a Slater determinant leading 
to the Fermi-Dirac distribution. 

In conclusion we have established a fundamental re- 
lation between adiabatic geometric phase and response 
function for cyclic as well as non-cyclic evolutions of 
chaotic quantum systems where chaos appears as a re- 
sult of the mean-field seen by a single-particle due to 
other particles. Since the imaginary part of the suscepti- 
bility is related to energy dissipation, our result provides 
a quantum mechanical way to interpret and understand 
decay of collective excitations of Fermi systems. The re- 
sults of this Letter can be generalised for the case when 
the dynamics is not fully chaotic using the results of [17]. 
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